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Abstract 



order phase transitions are initiated by the nucleation of bubbles of the stable phase. 

Traditionally, this process is described by the semiclassical nucleation theory devel- 

i-A . oped by Langer and, in the context of quantum field theory, by Callan and Coleman. 

G , They have shown that the nucleation rate T can be written in the form of the Ar- 

O . rhenius law: T = Ae~'^''. Here Tic is the energy of the critical bubble, and the 

prefactor A can be expressed in terms of the determinant of the operator of fluctu- 

^ ■ ations near the critical bubble state. It is not possible to find explicit expressions 

j^ ' for the constants A and Tic in the general case of a finite difference rj between the 

energies of the stable and metastable vacua. For small rj, the constant A can be 

determined within the leading approximation in r], which is an extension of the 

"thin wall approximation" . We have calculated the leading approximation of the 

prefactor for the case of a model with a real-valued order parameter field in two 

dimensions. 



*Poster presented at the "XIII International Congress on Mathematical Physics", July 17-22, 2000, 
London, UK. 



1 Introduction 

The problem of the decay of the metastable false vacuum at first order phase transitions 
has attracted considerable interest due to its numerous relations with condensed matter 
physics [0, quantum fields 0, cosmology |^, and black hole theory |^. In hanger's theory 
of homogeneous nucleation p , the false vacuum decay is associated with the spontaneous 
nucleation of a critical bubble of a stable phase in a metastable surrounding. In the 
context of quantum field theory, the nucleation theory was developed by Voloshin et al. 
0, and Callan and Coleman [0, H. The quantity of main interest is the nucleation rate 
r per time and volume. In the homogeneous nucleation theory it has the form of the 
Arrhenius law: 

r = ^exp(-K), (1) 

where Tic is the energy of the critical bubble. The prefactor A is determined by fiuc- 
tuations near the critical bubble state and can be expressed in terms of the functional 
determinant of the fiuctuation operator p[. 

It should be noted that in this article the so-called kinetic prefactor, which depends on 
the detailed non-equilibrium dynamics of the model, (see [|I|) is not included in F, which 
is therefore equal to twice the imaginary part of the free energy density of the metastable 
phase. 

In the general case, it is not possible to find the explicit critical bubble solution of 
the field equations analytically. However, the problem becomes asymptotically solvable, 
if the decaying metastable state is close enough in energy to the stable one, i.e. if the 
energy density difference rj between the metastable and stable vacua is small. The leading 
approximation in this small parameter is usually called the "thin wall approximation" 0, 
since at r/ ^ the critical bubble radius goes to infinity and becomes much larger than 
the thickness of the bubble wall. 

In the thin wall approximation, the critical bubble energy Tic can be easily obtained 
from hanger's nucleation theory. It turns out to be much more difficult to find explicitly 
the prefactor A in (|^). This problem, which is important for applications of nucleation 
theory, has been extensively studied in different models. 



A remarkable result on this subject was obtained by Voloshin |jTO|- He considered 
scalar field theory in 2 dimensions with a potential U{(f)) of the type shown in Fig. 1. 
Voloshin claimed that in the limit r/ ^ the nucleation rate F in such a model can be 
described by the simple universal formula 

F = ^exp(-K), (2) 

where the critical bubble energy Tic is given by 

«. = — . (3) 

V 

2 




Figure 1: The potential U with the false 



and true 



vacuum 



Here a is the surface tension of the wall between the stable and metastable vacua in the 
limit rj —>■ 0. Thus, according to |TD|, in this limit the nucleation rate F is determined 



by two well defined macroscopic parameters rj and a. Another claim of |10| is that there 
are no corrections to formula (0) and @ proportional to powers of the dimensionless 
parameter 77/cr^. Voloshin arrived at these conclusions by an analysis performed in the 
thin wall approximation. He replaced the original scalar field theory by an effective one, 
which describes only fluctuations of the critical bubble shape. This approach implies that 
all other fluctuations of the original scalar field could be properly accounted for by the 
correct choice of the macroscopic parameters 77 and a. 

Recently an analytical method was developed [0, which allows one to study nucle- 
ation in the scalar field model beyond the thin wall approximation. In |TT]] this method 
was used to calculate the nucleation rate for the first order phase transition in the three- 
dimensional Ginzburg-Landau model. In the present paper we apply the same approach 
to the two-dimensional case. We calculate the nucleation rate beyond the thin wall ap- 
proximation and verify directly Voloshin's claims (0) and @. 

Nucleation theory in two-dimensional scalar field theory has also been studied by 
Kiselev and Selivanov [|r^ , Strumia and Tetradis |TE[ , and other authors. In these articles, 
however, different renormalization schemes have been used and F has not been expressed 
in terms of macroscopic parameters 77 and a. This makes it difficult to compare their 
results with Voloshin's claims. 

In the article \li\ the nucleation rate was calculated in the two-dimensional Ising 
model in a small magnetic field. If Voloshin's results (|]) and (|^) are universal, they 
should be applicable as well to the Ising model in the critical region. Indeed, expressions 
(19), (23-26) of H^ rewritten in terms of t] and a are in a very good agreement with (^, 
d^). The exponent factors are the same, and the prefactors differ only by the number 
7r^/9 ~ 1.0966, which is very close to unity. This small discrepancy increased our interest 



in the subject of the present study. 



2 Model and notations 



We consider the two-dimensional asymmetric Ginzburg-Landau model defined by the 
Hamiltonian: 



7^(0) = £x 



l(9,0(x))2 + f/(0(x)) 



(4) 



where 0(x) is the continuous one-component order parameter, and the potential f/(0) 
depicted in Fig. 1 is given by 



2v 
Here Us{(j)) denotes the symmetric part of the potential: 



Us{<P) = f (0' - 



(5) 



(6) 



The potential U{(j)) has a metastable minimum (false vacuum) at (/) = </>+ and a stable one 
(true vacuum) at = 0_. The constant term Uo in (|^) is chosen to ensure f/(</>+) = 0. 
The partition function is given by the functional integral 



DcP exp [-n{^)] 



(7) 



The inverse temperature factor has been absorbed into 7i. 

It is convenient to define the mass m and the "inverse temperature" (3 = 1/T param- 
eters by 



d^ 



m 



■Usi^) 



gv 



P 



3m 



and to introduce dimensionless variables 



m 



""^ 2 ^'^' '^~~ 2m 



9 f~^ <P{x) 



V 



V p 



In dimensionless variables the Hamiltonian and partition function take the form 



l-i{^) = d'x 



-{Vvf + U{^{x)) 



where 



and 



m = 2 



(^^-l)'-K-l) 

Z = /d^(x) exp \-pn{ip) 



+ ^v{f-f+) 



(8) 



(9) 



(10) 



:iii 



(12) 



3 The critical bubble solution 

The uniform solutions of the field equation 

6H/6ip{x) = 0. 
are the stable (p- and false (metastable) ip+ vacua given by 



\^3 



^± 



3 6 27 ^ ' ^ 



(13) 



(14) 



The critical bubble v?fe(x) is the non-uniform radially symmetric solution of (0) approach- 
ing the false vacuum at infinity. That is, 



-p,^ - --7^ + 2<fh{<fb - 1) + o ^ 

ar'^ r ar 3 

Jim (pb{f^) 



0, 



(15) 



where f = Jx^x^. The profile of the critical bubble solution is shown schematically in 
Fig. 2. If fj is small, the thin wall centered sX f = R divides regions of false and stable 
vacua outside and inside the bubble, respectively. 




Figure 2: Profile of the critical bubble 



Equation (pIS]) can not be solved explicitly. Following the approach introduced by 
Miinster and Rotsch [jll| we shall construct the solution by expansion in powers of fj. 
Introducing the new independent variable ^: 



5 



(16) 



we write the following series for R and ^ptiO'- 



R 

MO 



V 



«o + aif] + a2ri + 0(7] ) 



M0 + M0v + M0r + o{v'). 



After substitution of (|T^-|TB|) into (^) one obtains perturbatively in rj: 



a-i 



-, ao = 0, ai 



ipoiO = tanh^, ipi{^) 



9' 
1 



(32 = 0, 



MO = - ^^ ^^^^2 I lOe - lee cosh(20 - 2e cosh(40 + 

2 In [2 cosh^] [^ + 8 sinh(20 + sinh(40] - 24 /rft ttanht 



(17) 
(18) 

(19) 



The bubble energy E = H [ipb{x)] can be written as 



E = 7r d^R + O 



-R 



f dMO X 



Substitution of (|T9]) into (|20|) yields 

E = 27i 



3?7^^U4 9 



oiv 



^3\ 



(20) 



(21) 



It is the basic principle of homogeneous nucleation theory that the decay of the metastable 
vacuum occurs through nucleation of the critical bubble. Callan and Coleman expressed 
the nucleation rate F of the metastable vacuum in terms of functional determinants [0, ^. 
In our notation their result takes the form 



/3E 1 



2-V^ 



.exp(^-f3E + S). 



(22) 



Here F = 4F/m^ is the dimensionless nucleation rate, and the entropy S associated with 
the critical bubble is given by 



expS" 



det;ww 

det'M 



1/2 



(23) 



where A^ and A^*^*^) are the fluctuation operators near the bubble fb{x) and metastable 
uniform vacuum </?+, respectively: 



M 



~m2 



-d' + 6 [Mr)] 
-.9^ + 6(^^-2. 
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(24) 
(25) 



The operator Ai has two zero modes proportional to d^(pb{x), fi = 1, 2, and one negative 
mode with the eigenvalue 

Ao = -4r7l (26) 

The notation det' implies that the three above mentioned modes are omitted in the 
corresponding determinant. After substitution of ([21| ) and (|^), equation (^2]) simplifies 
to 

In the subsequent sections we shall calculate the small fj expansion for the critical bubble 



entropy (^31) with accuracy o{rf). 



4 The bubble entropy 

Let us expand the bubble entropy S into the sum over the angular quantum number /i: 

oo 

S= Y. S,. (28) 



/x=— oo 



where 



|E^(lnAfj-lnA,-^) for /. ^ 0, ±1 

'~^ |lnAi°j + |g (inAfj-lnA,,) for/x = 0,±l. 

Here Aj^ and A^-^ are the eigenvalues of the radial Schrodinger operators corresponding 
to M and A<(°): 



(P' Id fi 



2 



^(0) ^ __^ ^ tL^Q^p'^ _2 (31) 

H' ^ ly Zi ly ri ly ty Zi ' 



Id u^ 

___ + ^ 

The eigenf unctions ^j^, Xjm ^-^e supposed to be finite at the origin f = and to obey the 
zero boundary condition at some large L : ipj^{L) = Xjfj.{L) = 0. 

4.1 Calculation of the "partial" entropy S^ 

Let us first consider the case ^ ^ 0, ±1. For such a \i the operators H^ and if^°^ are 
positive. Therefore the "partial" entropy Sp, can be written in terms of the trace of the 
resolvent operators: 

Ci 



where 



A{\, fi) = Sp 



(A -i/.)-^- (A -<))-' 



The integration path Ci shown in Fig. 3 hes in the right complex half-plane of A and goes 
around the spectra of both operators H^ and H^^ in the clockwise direction. It should 
be noted that the representation ( p2D is valid also for /i = 0, ±1. In the latter case the 
lowest discrete mode Ao^^ < becomes negative (/x = 0) or zero (/x = ±1). These modes 
contribute neither to (^) nor to (^). 

After deformation of the path Ci into the path C2 which goes around the logarithmic 
cut positioned at the negative real axis (see Fig. 3), one obtains 



1 r 

S, = -- JdXA{X,fi). 

—00 

In this equation one can proceed to the thermodynamic limit L ^ 00. 



(33) 




Figure 3: Integration paths Ci and C2 in the complex A-plane 



We obtained an exact representation for the integrand in (p3|). To describe it some 
notations are necessary. 

Let fi{r, A), gi{r, X), z = 1, 2 be the solutions of the linear ordinary differential equa- 
tion 

H^ ^(f) = A V'(r) (34) 

determined by their asymptotics: 



/i(r, A) -^ K^{qf), and /i(f, A) -^ /^(gr), at r ^ 00, 
5fi(f, A) -^ /C^(g_f), and5fi(f, A) ->/^(g_f), at f ^0. 



(35) 
(36) 



Here K^{z) and /^(-z) are the modified Bessel functions, the parameters q and g_ are 
defined as: 



q = (6y.^ - A - 2) , g_ = (6(^2. - A - 2) 



(37) 



Since the second order equation ([3^ ) has two hnearly independent solutions, there is a 
hnear dependence between the functions gi{r, A) and /j(f, A): 



3=1,2 



(3^ 



The above relation defines the 2 x 2-matrix a(A). The integrand in (|33|) can be expressed 
explicitly in terms of this matrix: 



A(A,/x) 



3^((/?2 



V- 



(A + 2-6v?i)(A + 2-6(^^ 



+ 



(ilna22(A) 
rfA ' 



(39) 



where A is assumed to be real and negative. 

The representation (|39D for the trace of resolvent operators is exact. However, equa- 
tion (p^ ) can not be solved in closed form for arbitrary fj. So we have to consider the 
small-r; expansion for 0:22 (A). We have obtained two terms of this expansion by use of 
a perturbation theoretical analysis of the scattering problem (|3^ - [38| ). Omitting the 
details, the logarithmic derivative of the matrix element 022 (A) up to quadratic terms in 
fj takes the form: 



din 0:22 (A) 
dX 



(4 + p2_A)i/2 



2p2 



+ 



{X-Ay A -4 A-3-p2 X-p2 



+ 0{fi'). (40) 



Here p is the angular momentum parameter defined asp = 2r]fi ^ Rfi, which is considered 
to be a quantity of order 1 for later purposes. 
Substitution of (|g) and (^ into (||) yields: 

2p2 



^(A,/i) 



1 



{4 + p^-Xf^ 

2\p\ 



{X-AY 



+ 



1 



+ 



A-4 A-3-p2 A-p- 



Integrating 



(A -4)^ 
we obtain the "partial" entropy S"^: 



(41) 



g ^H_(p' + 4) 



1/2 



4 



4 



1, 
+ -In 
2 



,1/2 



{p^ + AY'' + 2 (p2 + 4)1/2,^1 



p/ 



(p2+ 4)1/2 



(42) 



4.2 Summation of the /i-series 



Expression (^4 ) has a logarithmic singularity at p = 0. This singularity arises from the 
lowest discrete level of the Hamiltonian (0). To account for the contribution of this level 
properly, let us rewrite (pH]) as 



S= Y. s^{z)+\nR{z) 



(43) 



/i=— 00 



where 



s^^{z) 



\p\ (p^ + ^f 1, f[(p'+ 4)1/2 + 2)1' ^^2 + 4)1/2 + 1 



p2 _|_ ~/l 

1 °° n^ I 2 

^ /i=-oo /^ 



(p2+ 4)1/2 _1 



(44) 
(45) 



and z is some fixed positive number. 

First, we shall determine Ri^z). Taking into account the presence of the negative and 
zero modes, the following corrections are necessary: 



lni?(z) 



.2 , .2 1 ~ 477V + ^' 

o L In 



1 ^ p^ + ;2 

- > In ?; 

fi=— 00 f 



2^^^'" 4r/V 



1 ^ , 4f/V' + ^' 1, 2 W 2 . ~2N 



(46) 



These corrections arise from two facts: 

1. The logarithmic singularity of (|4^) at p ^ comes from the first discrete level 
Ao/i = 4f/^(/i^ — 1) of the operator (|^). 

2. The three eigenvalues Aq = — 47^^, A±i = of the operator (pOf ) already accounted 
for in (pTl) should be omitted in 5* and lni?(z). 

After these corrections we obtain from (|46| ) 



00 

R{z) =z(f+ Af) n 



/i2 + {z/2f,f 



where 



/.=2 /^ 

00 ,.2 I „2 



(2^)^^ 1 + ^ /(4)' 



/(^) = n 



2f/ \ 4772 
/i^ + x^ 2 sinh (ttx) 



2Ti' 



^^2 /^" - 1 TT x(x2 + 1) 



Therefore, 



(47) 



fl(.) = ^(2^)Mnh(.|.)^i(2,)fexpg), 
since fj is small, and z is a fixed number. Thus, we obtain 

In R{z) = 3 In (2r/) - In vr + — . 

This result can also be obtained by the ^-function technique described in [11]. 

For the calculation of the first sum in the right-hand side of ( ^3] ) it is convenient to 
apply Poisson's summation formula: 

00 00 °? 

E Sf,{z)= Yl I dfi s^{z)exp{2mfik). (48) 



k=-oo _oo 
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The integrals in terms with k ^ on the right-hand side converge. One can easily show 
that these integrals (at k ^ 0) are of order fj and can therefore be neglected in the adopted 
approximation. So, it is sufficient to calculate the k = term only, which we denote by 
s{z). Straightforward calculations yield: 



s{z 






(49) 



where S^ is given by (^). 

Now we can substitute (|49|) and ( ^7|) into (|43|) and obtain 



~\3 



exp^ = i^exp5(-i) 



TT 



1 + o(?7") 



(50) 



where 



5(-i)= Jdfis.^r'- 



As p — ^ oo, S*^ ~ — (2]9)^^. Therefore the integral in /i on the right-hand side diverges loga- 
rithmically at large |/i|. Indeed, this is the well-known ultraviolet divergency. Introducing 
a cut-off Pm in \p\ one finds: 



J 2f^ ^ 2f^\ ^"^ 2 V3, 
Relations (|50D , (|5lD are our final results for the bubble entropy. 



(51) 



5 Decay rate 

Substitution of ( pO] ) into (|27| ) yields for the dimensionless decay rate F: 

f = ^exp(-/?E(r/,r)), 



where 



37r 



£;(77,T) = |!^ - T A5(-^) + 0(77°) 

1 /27r T, , 11 TT , 



;^0n 



(52) 



(53) 



o{v 



To compare (0), (|53|) with Voloshin's result (H), (^), let us rewrite the latter in the 
dimensionless variables @ for the model (p!0|)-(|T^): 



Btt 



exp 



37ra^/5 

MM 



(54) 
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where a is the surface tension, and (</?) is the average value of the scalar order parameter 
at r] = 0. In (0) we have replaced the energy difference between the false and true vacua 
by I fi {if) (3. 

It should be stressed that the parameters {(p) and a in (|5^) are exact macroscopic 
quantities which depend on the interaction parameter g, or in other terms, on the "tem- 
perature" T = [3~^ = g/{?>rn?). At low "temperatures" these quantities can be expanded 
in T: 

(^) = 1 + Tp + 0(T2), a=^+Td + 0{T^), (55) 

and the corrections linear in T can be found perturbatively: 

3/1 \ 

^ = _31np^_3^ 1 

Here i) is the one-loop correction in (</?), the correction in a was calculated in the strip 
geometry from the ratio of partition functions j3a = L^^ \yi[Zq/Zi) with different boundary 
conditions: 

Zi : -V9(xi +Li,X2) = V5(xi,X2 + L2) = V5(xi,X2). 

It should be noted that both p and d contain logarithmic ultraviolet divergencies and 
depend on the cut-off procedure. In (|56D , (|57D we have chosen a similar cut-off-procedure 
as in d^). Substituting (p5D-(p7|) into ( pH ) and expanding in T one obtains 



4/5/7 
^— exp 



2vr 1 / , 11 TT ^, ' 

— 51np™ + — --= + 0(T) 



3Tf/ 277 V 2 V3 



in agreement with our results (|52D , (|53D . 



Summary: Our semiclassical calculation of the nucleation rate T in the two- 
dimensional Landau- Ginzburg 0^-model confirms Voloshin's results (|^), (^), which were 
derived in the thin wall approximation. In particular, we confirm the prefactor value 
A = ri/{27i) first obtained by Kiselev and Selivanov |12|, and Voloshin [[l^ . 

This value differs from that obtained for the two dimensional critical Ising model M 



by the numerical factor 7r^/9 ~ 1.0966. We suppose that this small discrepancy is the 



result of approximations used in |JT^, and the prefactor value A = ?7/(27r) is universal. 
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